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Abstract 

The seesaw mechanism in models with extra dimensions is shown to be generically consistent 
with a broad range of Majorana masses. The resulting democracy of scales implies that the seesaw 
mechanism can naturally explain the smallness of neutrino masses for an arbitrarily small right- 
handed neutrino mass. If the scales of the seesaw parameters are split, with two right-handed 
neutrinos at a high scale and one at a keV scale, one can explain the matter-antimatter asymmetry 
of the universe, as well as dark matter. The dark matter candidate, a sterile right-handed neutrino 
with mass of several keV, can account for the observed pulsar velocities and for the recent data from 
Chandra X-ray Observatory, which suggest the existence of a 5 keV sterile right-handed neutrino. 
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I. INTRODUCTION 



The possible existence of extra dimensions has changed the way one thinks about the 
hierarchy of scales in particle physics If the light Standard Model (SM) degrees 

of freedom are localized on a (3+l)-dimensional brane embedded in a higher-dimensional 
spacetime, while any other degrees of freedom reside on other branes or propagate in the 
bulk, the perceived difference between the high scale and the low scale may arise from the 
suppression that is exponential in the size of the extra dimension. It is a generic phenomenon 
in this class of models that the interactions of fields localized on different branes can be 
exponentially suppressed by the distance between the branes j^] . In the low-energy effective 
theory, the smallness of the couplings may seem completely unnatural, because the symmetry 
group of the (3+l)-dimensional effective theory need not increase in the limit where the size 
of the extra dimension becomes large. 

Neutrino masses and the possibility of leptogenesis offer an intriguing connection to the 
high-scale physics. In (3+1) dimensions, the naturalness arguments suggest that the SU(2) 
singlet fermions should have masses of the order of the high scale, in contrast with the non- 
singlet fermions, whose masses are protected by the chiral gauge SU(2) symmetry. However, 
if extra dimensions exist, one must rethink the naturalness arguments. A fermion localized 
on a distant brane may appear as weakly interacting light particle in the low-energy effective 
four- dimensional Lagrangian. 

Let us review the natural values of parameters in the modern version of the Standard 
Model, by which we mean the original Standard Model of Glashow, Salam, and Weinberg, 
supplied with three right-handed neutrinos generating the observed neutrino masses via the 
seesaw mechanism {4]. As mentioned above, and as we will show in detail below, both heavy 
and light right-handed neutrinos can be equally natural if extra dimensions exist, assuming 
that the small-mass fermions in (3+1) dimensions arise from the heavy particles on a remote 
brane. The simultaneous existence of light and heavy right-handed neutrinos allows one to 
explain both the baryon asymmetry of the universe and dark matter in the framework of 
the simplest model consistent with the data, - the modern variant of the Standard Model. 
The dark matter candidate in this model is a keV sterile neutrino, the possible existence of 
which is supported by some astrophysical data 5|-[8|]. 

The only way to generate the baryon asymmetry of the universe in this minimal scenario 
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is via leptogenesis js], because the lack of CP violation and the weakness of the electroweak 
phase transition eliminate any alternative scenarios. Only two right-handed neutrinos are 
necessary for leptogenesis [lo| . However, having three right-handed neutrinos, one for each 
generation of fermions, is arguably more appealing; this also allows one to embed the Stan- 
dard Model into SO(10) Grand Unified Theory or any theory that has a gauge U(1)b-l 
symmetry. What is the natural value for the Major ana mass of the third right-handed field? 
In the presence of extra dimensions, the high-scale values are not necessarily favored by 



any naturalness arguments. However, if this Majorana mass is of the order o 
can explain dark matter, as well as the long-standing puzzle of the pulsar kicks 



same particles can play an important role in the formation of the first stars [3] and other 



ceV, one 
The 



astrophysical phenomena ll|. Furthermore, at least two recent astrophysical observations 
suggest the existence of a sterile neutrino either with 5 keV mass js| or 17 keV mass 12]. 

By splitting the scales of the heavy and light right-handed neutrinos, which is naturally 
realized in the split seesaw model described below, we achieve an elegant and simple de- 
scription of all known experimental data. We note that models with three right-handed 

n n 

neutrinos below the electroweak scale have been proposed [13|, |14J, and in one of them, 



dubbed z/MSM 



14j . one can explain the baryon asymmetry by demanding a high degree 



of degeneracy between two Majorana masses, both of the order of several GeV. The high 
degree of degeneracy is required to amplify the effects of CP violation in a leptogenesis sce- 
nario that involves oscillations 15]. In contrast, we employ the standard leptogenesis using 
decays of the heavy right-handed neutrinos, and no mass degeneracy is required. We em- 
phasize, however, that our model with an extra dimension is not limited to the conventional 
leptogenesis scenario, and it can explain a very broad range of scales. Indeed, our model 
assures that the beauty of the seesaw formula, which explains the smallness of the neutrino 
masses by relating them to the ratio of the weak scale and the GUT scale, is preserved for 
practically arbitrary choice of the right-handed neutrino masses. 
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II. SPLIT SEESAW MECHANISM 



We begin by reviewing the seesaw mechanism in the 4D theory. The relevant terms in 
the Lagrangian are given by 

C = iNtfdpNi + (\iaNiLa <t> - ^M ffi iVfiV 4 + h.c.) , (1) 

where TVj, L a and <fi are the right-handed neutrino, lepton doublet and Higgs boson, respec- 
tively, i denotes the generation of the right-handed neutrino, and a runs over the lepton 
flavor, e, /i and r. Integrating out the massive right-handed neutrinos yields the seesaw 
formula for the light neutrino mass: 

(™»U = E A ^^" ( 2 ) 



The atmospheric 



16| and solar 
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neutrino oscillation experiments have provided firm 
evidence that at least two neutrinos have small but non-zero masses, and the mass splittings 
are given by Am 2 tm ~ 2 x l(T 3 eV 2 and Am| ~ 8 x l(T 5 eV 2 . The seesaw mechanism then 
suggests that a typical mass scale of the right-handed neutrinos should be ~ 10 15 GeV, close 
to the GUT scale, for Aj Q ~ 1. 

To explain the neutrino oscillation data, one must introduce more than one right-handed 
neutrino. It was shown in Ref. lOj that two right-handed neutrinos, N 2 and N 3 , would suffice 
for this purpose. Moreover, it is possible to generate the cosmological baryon asymmetry via 
leptogenesis with the two right-handed neutrinos. Thus, the addition of the two right-handed 
neutrinos appears to be the most economical extension of the SM. 

However there are two issues in the above model, which, as we will see below, can be 
addressed simultaneously and naturally in a theory with an extra dimension. One issue 
is that at least the lightest of N 2 ^ must have a mass OflO 11-12 ) GeV, several orders of 



magnitude below the GUT scale, for the leptogenesis to work [19l l20l|. Therefore, to explain 
the neutrino masses, one must suppress both the right-handed Majorana mass and the 
associated Yukawa couplings. This is not a severe fine-tuning, and it can be rectified by the 
Froggatt-Nielsen mechanism {2]] , but the introduction of a new symmetry and the breaking 
of this symmetry at an appropriate scale, for this purpose, appear somewhat ad hoc. The 
other issue is that, if one gauges the U(1)b-l symmetry, three right-handed neutrinos are 
required to achieve the anomaly cancellation. Then the question arises: what is the role of 



the third right-handed neutrino, iVi? In one limiting case, N x may have a Planck-scale mass, 
and this particle would play no role in the low energy physics. In particular, it is clear from 
the seesaw formula that such a heavy N% does not contribute to the light neutrino mass, m v . 
In the other limit, N\ may be very light. An interesting possibility arises in this case; if the 
mass is at a keV scale, then N\, which is long-lived on cosmological time scales, is a viable 
dark-matter candidate 22| . In addition, the same particle would be produced in a supernova 



explosion, and it would be emitted with an anisotropy of a few per cent, which is sufficient 



to explain the observed velocities of pulsars [6|. However, one could ask whether such an 
extremely light mass may upset the seesaw mechanism, and it is also unclear how such a 
small mass may arise naturally. We will show that, in a theory with an extra dimension, a 
split mass spectrum arises naturally, without disrupting the seesaw mechanism. The split 
seesaw mechanism can explain both mild and large mass hierarchy, i.e., why < Mgut 
and why M x < M 2j3 . 

Let us consider a 5D theory compactified on S 1 /Z 2 with coordinate y G [0, £}. One of the 
boundaries at y = is identified with the SM brane, where the SM degrees freedom reside, 
while the other at y = I is a hidden brane. The size of the extra dimension i and the 5D 
fundamental scale M are related to the 4D reduced Planck scale as 

Ml = M 3 £. (3) 

For the consistency of the theory, the compactification scale M c = \ jl must be smaller than 
M. As reference values we take M ~ 5 x 10 17 GeV, M c ~ 10 16 GeV and the B — L breaking 
scale Vb-l ~ 10 15 GeV in the following. 

We introduce a Dirac spinor field in 5D, \1/ = (x a ,'0 a ) T , with a bulk mass m: 

S = J d A xdyM (i^r A d A ^ + mM) , (4) 

where A runs over 0, 1, 2, 3, 5, and the 5D gamma matrices T A are defined by 

(5) 

We have inserted the mass scale M in Eq. (j3J) so that the mass dimension of \I/ is 3/2 as in 
the 4D case. The zero mode of \& should satisfy 

(iT 5 d 5 + m)¥ 0) = 0. (6) 




The bulk profile of the zero mode is therefore given by exp(=pm?/) for \ and tjj [23|, [24 1. 

A 4D chiral fermion can be obtained from the 5D Dirac fermion by orbifolding. If we 
assign a Z 2 parity —1 and +1 to x and ip, respectively, we can see that only ip has a zero 
mode with an exponential profile in the bulk. (For consistency we have assigned a negative 
Z 2 parity to the bulk mass m.) The zero mode of ty R = (0, ip) T can be expressed in terms 
of the canonically normalized (right-handed) fermion ip R D ^ in the 4D theory as 



•8W) = ^^^e^r> W . (7) 

As one can see from the ^-dependence of the wave function, the zero mode peaks at the 
hidden brane at y = £ for positive m. In particular, in the limit of mi 3> 1, the overlap of 
the zero mode with particles on the SM brane becomes exponentially suppressed. 

We would like to promote the right-handed neutrino A 7 , to a bulk field in 5D 3J. For 
definiteness, we identify the zero modes of with the right-handed neutrinos in the 4D 
theory. Then the zero modes have an exponential profile in the y direction ~ exp(my). We 
also introduce a U(1)b-l gauge field in the bulk and a scalar field $ with a B — L charge 
—2 on the SM brane. Assuming that the $ develops a VEV ^b-l ~ 10 15 GeV, the zero 
mode of the U(1)b-l gauge boson receives mass via the usual Higgs mechanism, and there 
is no additional zero mode. The VEV of $ also gives rise to heavy Majorana masses for the 
right-handed neutrinos. After integrating out the heavy $ and the U(1)b-l gauge boson, 
one obtains the Lagrangian for the zero modes: 

S = J d*xdy {M (itfgr^g + mM°Mi) 

+%) (^B-L*SX(°) + X^fX + h.C.) } , (8) 

where Ki and Aj Q are numerical constants of order unity, and we introduced the lepton and 
Higgs doublets on the SM brane at y — 0.* 1 Using the relation (171), we obtain the effective 



4D mass and Yukawa couplings in ([TJ as 26] 



, r 2rrii 

M Ri = fi^B-L 777 om-f TT' ( 9 ) 



'M{e 2m > e - 1 



\ia I 2m*, r / M Ri 




(10) 



* 1 Since all the right-handed neutrinos are in the bulk, the U(1)b-l becomes anomalous on the SM brane. 
However, one can cancel the anomaly by introducing Chcrn-Simons terms with an appropriate coefficient 
in the bulk [25|. 
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As expected, the right-handed neutrino mass and the Yukawa couplings are suppressed by 
an exponential factor. In the extreme limit of ^> 1, both Mr and A are extremely 
suppressed. What is remarkable is that the exponential factors cancel in the the seesaw 
formula for the light neutrino masses ([2]) because the Yukawa coupling squared is suppressed 
by the same factor as the right-handed mass Thus, the neutrino masses are given by the 
familiar seesaw relation: 




where the quantity in the parenthesis is of order unity. We can clearly see that no small 
parameters such as exponentially suppressed mass or coupling appear in the seesaw formula; 
the typical neutrino mass scale is given by the ratio of the square of the weak scale to 
the B — L breaking scale, showing that the seesaw mechanism is robust against splitting 
the right-handed neutrino mass spectrum. So we can easily realize a split mass spectrum 
by choosing appropriate (not extremely large or small) values of rrii, without spoiling the 
seesaw mechanism. For instance, we obtain Mr 2 ~ 10 12 (10 11 ) GeV and Mjh ~keV for 
ni2 — 2.3(3.6)^ _1 and for ni\ ~ 24£ _1 , respectively, where we set Ki — 1 and used the 
reference values of v-b-l, £ arid M.* 2 



III. COSMOLOGY 



In the context of sterile neutrino dark matter, it is customary to parametrize it in terms 
of the effective mixing angle 9 and the mass m s . They are given by 

m s = M RU (12) 
K 3 x 10 -» | '^EJVI^ ( 10^GeV\ (M\ (13) 



m u ~ 2 m s , (14) 

where we have assumed 9 <C 1. If Ni is the dark- matter particle, the X-ray constraints 
imply 9 2 < 10~ 5 (m s /keV) -4 , while the small-scale structure constraints imply m s > 2keV 



# 2 It is also possible to realize Mr^ — Mm ~ GeV for ni2 — ~ 171 1 , as considered in the i/MSM 14|. 
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for the resonant production [27| (see, e.g., Refs. [5|, |28[ for review). For the above range 



of parameters the contribution of JVi to the light neutrino mass is negligible 



29|. The 



prospects of direct detection of sterile dark matter are also stymied by the smallness of 



the mixing angle 



30] . In contrast, the X-ray instruments can search for decays of sterile 



right-handed neutrinos, which have the half-life much longer than the age of the universe. 
Here the smallness of sin 2 9 is compensated by the huge number of dark-matter particles. In 
particular, if the 2.5 keV spectral feature in the recent Chandra X-ray Observatory data 8| is 
explained by the decay of a 5 keV relic sterile right-handed neutrino, the inferred parameters 
m s ~ 5keV and 9 = (0.2 — 1.4) x 10~ 9 js] imply that such a sterile neutrino does not 
contribute to the neutrino mass. 

To account for the observed dark matter density, sterile right-handed neutrinos must 
have the correct abundance. A generic way to produce relic sterile neutrinos is through non- 
resonant (NR) oscillations 22]. For the mixing angle sin 2 9 ~ 10~ 9 and the mass m s ps 5keV, 



the resulting abundance of sterile neutrinos is in the right range |22l. l3lM33|. Assuming that 
the sterile neutrinos produced by the NR oscillations account for the total dark matter 
density, there is an upper limit m s < 5keV from the X-ray observations 34|. On the other 



hand, the lower limits m s > 8 keV 



35] and m s > 11 keV 36] based on the Ly-a observations 



and on dwarf spheroidals, respectively, appear to rule out such dark matter produced by NR 
oscillations. Thus, if the X-ray line observed by the Chandra X-ray Observatory is due to the 
decay of dark matter sterile neutrinos, they must be produced by some other mechanism. 
Such mechanisms, producing colder sterile neutrinos, consistent with all the constraints, 
have been studied in the context of other models [37]. We will show that, in the split seesaw 
scenario, the sterile neutrinos can be produced with the same or greater cool-down factor as 



those produced at the electroweak scale 



371 ] . Thus, our dark-matter candidate presents no 



conflict with the Lyman-a and other small-scale structure bounds. 

In our set-up, the sterile neutrino has a U(1)b-l gauge interaction, which may become 
important in the early Universe. In particular, the reheat temperature Tr must be greater 
than Mr2 or Mr 3 for the leptogenesis to work. According to Ref. [20[, a successful lepto- 
genesis is possible for Tr > 10 11 GeV, unless Mr 2 and Mr 3 are extremely degenerate. For 
such a high reheat temperature, the main production process is pair production of N\ from 
the SM fermions in plasma through the s-channel exchange of the B-L gauge boson. The 
production is most efficient at reheating. The number to entropy ratio of the sterile neutrino 



S 



produced by this mechanism can be roughly estimated as 

n Nl {av)n)/H 



Nl 
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T=T R 

3 



9* \ 2 / *>b-l \- 4 ( Tj 



ID" 4 TT^TT \ K „," r , (15) 



10 2 / Vl0 15 GeV/ V 5xl ° 13GeV 
where H is the Hubble parameter, g* counts the relativistic degrees of freedom at the 
reheating, (av) ~ T 2 /v^_ L is the production cross section, nf ~ T 3 is the number density 
of the SM fermions in plasma, and the first equality is evaluated at the reheating. The 
numerical solution of the Boltzmann equation gives a consistent result 
abundance is related to the mass density to entropy ratio as 



. The dark matter 



CW = 0.2 x (^-] ( Ym - ) . (16) 
dark V5keV/ V0.7 x 10" 4 J V ; 

Thus, the reheat temperature as large as 10 13 GeV is needed to account for all the dark 
matter density by this production process. The thermal leptogenesis works with such a 
high temperature. Also since the sterile right-handed neutrinos are out of equilibrium since 
their production at high temperature, when all the Standard model degrees of freedom are 
thermally excited, the average momentum of such dark-matter particles is red-shifted by a 
factor of a cubic root of the ratio of the degrees of freedom, which is about 3.5 (37 ]. 

Another scenario for producing dark matter with the correct abundance (but with differ- 
ent kinetic properties) places no constraint on the reheat temperature and predicts an even 
colder dark matter, as long as reheating restores the gauge U(1)b-l symmetry. We assume 
that the U(1)b-l symmetry is broken spontaneously at 10 15 GeV by the VEV of the Higgs 
boson $. Since the right-handed neutrinos are coupled to the {7(1)b-l gauge boson with 
the universal gauge coupling go, they are produced at high temperature and they reach the 
equilibrium density n eq (T) with a distribution function f cq (E, T) = l/(exp(E/T) + 1) above 
the transition temperature T t . However, after the phase transition, the U(1)b-l gauge boson 
becomes massive and the production of iVj is suppressed by the gauge boson mass (just as 
in the production mechanism discussed above). Meanwhile, the rest of the Standard Model 
particles are produced in the true vacuum. The entropy produced in this phase transition 
can dilute the number density of dark-matter sterile neutrinos from its equilibrium value 
Y^i = 45£(3) 2/ (2ir A g*) ~ 10~ 2 to what is required to explain dark matter, see eq. (Tl6|) . In 
this case, the reheat temperature is not constrained from above, but the final temperature 



after the phase transition Tf must satisfy the same constraint as the reheat temperature in 
the previous scenario. 

If the Higgs potential at zero temperature is 

^0) = -^ + ^, (17) 

where we fix the vev of the Higgs to be i>b-l> i-e., 

JL = Ub „ l w io 15 GeV. (18) 
v A 

The temperature-dependent effective potential [39j is 

T) « V(* 0) + ?^TV - 3g3 + 9 2 A 4 ; 3A3/2 ^. (19) 

where g = 2 c/ is the effective gauge (B — L) coupling of the Higgs boson, which has (B — L) = 
—2. The phase transition from the U(l) symmetric vacuum to the broken-symmetry vacuum 
takes place when the tunneling rate 39j, |40[ per Hubble volume per Hubble time is of order 
one: 

T/(^)V^ x H t -*~1, (20) 

where 

H (fgA 1/2 Tl (21) 

1 \ 90 ) M P 1 ; 



Here S3 is the action o 
approximately given by 



the three-dimensional bounce, and in the thick-wall limit, it is 



39] 

S 3 44( 7 (T 2 -T 2 )) 3/2 



T a 2 T 3 

with 



(22) 



_ 3g 3 + g\ + 3X 3 / 2 _ 3g 2 + AX jx_ 
a ~ 8tt ' 7 " 12 ' c " V7 



The correct abundance of dark matter is attained if the ratio of entropies before and after 
the phase transition is given by 

j|~^~10 2 , or I>~5T t . (23) 

To avoid overproduction of sterile neutrino after the phase transition, the final temperature 
should be smaller than the reheat temperature estimated in the previous scenario, 

T f < T max = 5 x 10 13 GeV (^f^) ^ • (24) 
10 



By energy conservation, 

T f ~ ( e /#*) 1/4 ~ 0.3e 1/4 , 

where e = — V(i>b-L) 0). Note here that the latent heat comes from the energy difference 
between the symmetric phase = and the broken phase = ^b-Lj which is estimated by 
using the zero-temperature potential. 

While there is a significant freedom in the choice of parameters, one can verify that, for 
example, the choice of g = 2go = 0.3, fi = 10 11 GeV satisfies the constraints (I23|24p for 
T t ~ 7 x 10 n GeV and T f « 5T t ~ 3 x 10 12 GeV. For an entropy dilution of 0(100) to 
occur, a relatively light mass [i is required, which leads to some degree of fine-tuning of the 
parameters in the scalar potential. This fine-tuning may be ameliorated by extending the 
Higgs sector, so as to increase the strength of the first-order phase transition. 

While the X-ray constraints depend only on the mass, mixing, and abundance of sterile 



neutrinos, the clustering properties of dark matter, probed by Lyman-a bounds [35| and the 
observations of dwarf spheroidal galaxies 36|,|4l|, are sensitive to the momentum distribution 
of dark-matter particles. In our first scenario, the limited preheating produces a non-thermal 
distribution of dark matter particles, which is further cooled down and redshifted in the 
subsequent history of the universe. In our second scenario, in addition to this redshifting, 
the B— L breaking phase transition redshifts the momenta of dark-matter particles by an 
additional factor Tf/T t m 5. Thus, the two scenarios produce dark matter with different 
free-streaming lengths, each of which is small enough to satisfy the present bounds. 



IV. DISCUSSION 



We have shown that the seesaw formula is robust with respect to splitting the right- 
handed neutrino mass spectrum in a model with an extra dimension. The Yukawa couplings 
and the Majorana masses are suppressed in such a way that the usual seesaw relation 
is preserved. This remarkable feature allows the use of the seesaw mechanism for two 
heavy right-handed neutrinos and one keV sterile right-handed neutrino in a model that 
does not require any unnatural fine-tuning of parameters. The resulting split seesaw model 
explains the origin of the baryon asymmetry of the universe, dark matter, and the smallness 
of the active neutrino masses in an natural and holistic manner. Furthermore, the light 
sterile neutrino provides a simple explanation for the observed velocities of pulsars, due 
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to anisotropic emission of such particles from a supernova. Finally, the model offers an 
explanation as to why there are three generations of fermions in nature: two right-handed 
neutrinos are needed for leptogenesis, and the third one plays the role of dark matter. X-ray 
astronomy offers an exciting opportunity to discover this dark matter candidate, and indeed 
the recent data from Chandra X-ray Observatory show a spectral feature consistent with 
this form of dark matter. If this result is confirmed, the detailed analysis of the structure 
on sub-galactic scales can help distinguish between the two production mechanisms we have 
discussed. 

Acknowledgments 

F.T. thanks Antonio Riotto and CERN Theory Group for the warm hospitality while 
part of this work was completed. A.K. was supported by DOE Grant DE-FG03-91ER40662 
and NASA ATP Grant NNX08AL48G. The work of F.T. was supported by the Grant-in- Aid 
for Scientific Research on Innovative Areas (No. 21111006) and Scientific Research (A) (No. 
22244030), and JSPS Grant-in-Aid for Young Scientists (B) (No. 21740160). This work 
was supported by World Premier International Center Initiative (WPI Program), MEXT, 
Japan. A.K. thanks Aspen Center for Physics for hospitality. 



[1] N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. Lett. B 429, 263 (1998) 
[arXiv:hep-ph/9803315] . 

[2] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 (1999) [arXiv:hep-ph/990522 1 1 ; L. Ran- 
dall and R. Sundrum, Phys. Rev. Lett. 83, 4690 (1999) [arXiv:hep-th/9906064| . 

[3] A. Hebecker, J. March-Russell and T. Yanagida, Phys. Lett. B 552, 229 (2003) 
[arXiv:hep-pfi/0208249] . 

[4] T. Yangida, in Proceedings of the "Workshop on the Unified Theory and the Baryon Number 
in the Universe", Tsukuba, Japan, Feb. 13-14, 1979, edited by O. Sawada and A. Sugamoto, 
KEK report KEK-79-18, p. 95, and "Horizontal Symmetry And Masses Of Neutrinos" , Prog. 
Theor. Phys. 64 (1980) 1103; M. Gell-Mann, P. Ramond and R. Slansky, in "Supergravity" 
(North-Holland, Amsterdam, 1979) eds. D. Z. Freedom and P. van Nieuwenhuizen, Print-80- 



12 



0576 (CERN); see also P. Minkowski, Phys. Lett. B 67, 421 (1977). 
[5] A. Kusenko, Phys. Rept. 481, 1 (2009) |arXiv:0906.2968l [hep-ph]]. 

[6] A. Kusenko and G. Segre, Phys. Lett. B 396, 197 (1997) |arXiv:hep-ph/9701311| ; 
G. M. Fuller, A. Kusenko, I. Mocioiu and S. Pascoli, Phys. Rev. D 68, 103002 (2003) 
|arXiv:astro-ph/0307267] ; A. Kusenko, B. P. Mandal and A. Mukherjee, Phys. Rev. D 77, 
123009 (2008) [arXiv: 080 1.47341 [astro-ph]]. 

[7] P. L. Biermann and A. Kusenko, Phys. Rev. Lett. 96, 091301 (2006); M. Mapelli, A. Ferrara 
and E. Pierpaoli, Mon. Not. Roy. Astron. Soc. 369, 1719 (2006); J. Stasielak, P. L. Bier- 
mann and A. Kusenko, Astrophys. J. 654, 290 (2007); Acta Phys. Polon. B 38, 3869 (2007); 
E. Ripamonti, M. Mapelli and A. Ferrara, Mon. Not. Roy. Astron. Soc. 375, 1399 (2007). 

[8] M. Loewenstein and A. Kusenko, Astrophys. J. 714, 652 (2010) |arXiv:0912 .05521 [astro- 
ph.HE]]. 

[9] M. Fukugita and T. Yanagida, Phys. Lett. B 174, 45 (1986). 
[10] P. H. Frampton, S. L. Glashow and T. Yanagida, Phys. Lett. B 548, 119 (2002) 
[arXiv:hep-ph/0208157] . 

[11] F. Munyaneza, PL. Biermann, Astron. and Astrophys., 436, 805 (2005); ibid, 458, L9 (2006). 
[12] D. A. Prokhorov and J. Silk. larXiv:1001.0215l [astro-ph.HE], 
[13] A. de Gouvea, Phys. Rev. D 72, 033005 (2005) [arXiv:hep-ph/0501039] . 
[14] T. Asaka, S. Blanchet and M. Shaposhnikov, Phys. Lett. B 631, 151 (2005) 
|arXiv:hep-ph/0503065 . 

[15] E. K. Akhmedov, V. A. Rubakov and A. Y. Smirnov, Phys. Rev. Lett. 81, 1359 (1998) 

[arXiv:hep-ph/9803255] ; T. Asaka and M. Shaposhnikov, Phys. Lett. B 620, 17 (2005). 
[16] J. Hosaka et al. [Super-Kamiokande Collaboration], Phys. Rev. D 74, 032002 (2006) 

arX iv :hep-ex / 06040 11] ; Y. Ashie et al. [Super-Kamiokande Collaboration], Phys. Rev. D 71, 

112005 (2005) [arXiv:hep-ex/0501064| . 
[17] J. Hosaka et al. [Super-Kamkiokande Collaboration], Phys. Rev. D 73, 112001 (2006) 

|arXiv:hep-ex/0508053| ; J. P. C ravens et al. [Super-Kamiokande Collaboration], Phys. Rev. D 

78, 032002 (2008) |arXiv:0803.43T2l [hep-ex]]. 
[18] B. Aharmim et al. [SNO Collaboration], Phys. Rev. C 72, 055502 (2005) 

[arXiv:nucl-ex/0502021] ; larXiv:0910.2984l [nucl-ex] . 
[19] T. Endoh, S. Kaneko, S. K. Kang, T. Morozumi and M. Tanimoto, Phys. Rev. Lett. 89, 



13 



231601 (2002) [arXiv:hep-ph/0209020| . 

M. Raidal and A. Strumia, Phys. Lett. B 553, 72 (2003) [arXiv:hep-ph/0210021] . 

C. D. Froggatt and H. B. Nielsen, Nucl. Phys. B 147, 277 (1979). 

S. Dodelson and L. M. Widrow, Phys. Rev. Lett. 72, 17 (1994) [arXiv:hep-ph/9303287] . 
R. Jackiw and C. Rebbi, Phys. Rev. D 13, 3398 (1976). 

D. B. Kaplan, Phys. Lett. B 288, 342 (1992) |arXiv:hep-lat /9206013| . 

K. I. Izawa, T. Watari and T. Yanagida, Phys. Lett. B 534, 93 (2002) |arXiv:hep-ph/0202171] . 



A. Hebecker and J. March-Russell, Nucl. Phys. B 613, 3 (2001) [arXiv:hep-ph/0106166] . 
A. Boyarsky, J. Lesgourgues, O. Ruchayskiy and M. Viel, Phys. Rev. Lett. 102, 201304 (2009) 
|arXiv:0812.3256l [hep-ph]]. 

A. Boyarsky, O. Ruchayskiy, M. Shaposhnikov, Ann. Rev. Nucl. Part. Sci. 59, 191-214 (2009). 
|arXiv:0901.0"0TT1 [hep-ph]]. 

A. Boyarsky, A. Neronov, O. Ruchayskiy and M. Shaposhnikov, JETP Lett. 83, 133 (2006) 



|arXiv:hep-ph/0601098] . 

S. Ando and A. Kusenko, Phys. Rev. D 81, 113006 (2010) |arXiv:1001.5273l [hep-ph]]. 
K. Abazajian, G. M. Fuller and M. Patel, Phys. Rev. D 64, 023501 (2001) 
|arXiv:astro-ph/0101524] ; A. D. Dolgov, Phys. Rept. 370, 333 (2002) |arXiv:hep-ph/0202122] . 
K. Abazajian, G. M. Fuller and W. H. Tucker, Astrophys. J. 562, 593 (2001) 



|arXiv:astro- ph/0106002| . 

T. Asaka, M. Laine and M. Shaposhnikov, JHEP 0701, 091 (2007) |arXiv:hep-ph/0612182] , 
M. Loewenstein, A. Kusenko and P. L. Biermann, Astrophys. J. 700, 426 (2009) 
|arXiv:0812.27T0l [astro-ph]]. 

A. Boyarsky, J. Lesgourgues, O. Ruchayskiy and M. Viel, JCAP 0905, 012 (2009) 



[ar Xiv:0812.0010l [astro-ph]]. 

E. Polisensky and M. Ricotti. [arXiv:1004.145"9l [astro-ph.CO], 

A. Kusenko, Phys. Rev. Lett. 97, 241301 (2006) [arXiv:hep-ph/0609081] . K. Petraki and 
A. Kusenko, Phys. Rev. D 77, 065014 (2008) jar Xiv:0711.4646l [hep-ph]]: D. Boyanovsky, Phys. 
Rev. D 78, 103505 (2008) |arXiv:0807.0646l [astro-ph]]: J. Wu, C. M. Ho and D. Boyanovsky, 
Phys. Rev. D 80, 103511 (2009) [ arXiv:0902.4278 [hep-ph]]. 
[38] S. Khalil and O. Seto, JCAP 0810, 024 (2008) jarXiv:0804.0336l [hep-ph]]; 
G. Gelmini, S. Palomares-Ruiz and S. Pascoli, Phys. Rev. Lett. 93, 081302 (2004) 



11 



arXiv:astro-ph/0403323 1; G. Gelmini, E. Osoba, S. Palomares-Ruiz and S. Pascoli, JCAP 

0810, 029 (2008) |arXiv:0803.2735l [astro-ph]]. 
[39] D. A. Kirzhnits and A. D. Linde, Annals Phys. 101, 195 (1976); A. D. Linde, Nucl. Phys. B 

216, 421 (1983) [Erratum-ibid. B 223, 544 (1983)]. 
[40] A. Kusenko, Phys. Lett. B 358, 47 (1995) [arXiv:hep-ph/9506386] ; 

A. Kusenko, K. M. Lee and E. J. Weinberg, Phys. Rev. D 55, 4903 (1997) 

|arXiv:hep-th/9609i00] . 

[41] G. Gilmore, M. I. Wilkinson, R. F. G. Wyse, J. T. Kleyna, A. Koch, N. W. Evans and 
E. K. Grebel, Astrophys. J. 663, 948 (2007) [arXiv:astro-ph/0703308| ; R. F. G. Wyse and 
G. Gilmore. laFXiv: 0708 .14921 [astro-ph] . G. Gilmore et aL larXiv:0804.1919l [astro-ph]: D. Boy- 
anovsky, H. J. de Vega and N. Sanchez, Phys. Rev. D 77, 043518 (2008); H. J. de Vega and 
N. G. Sanchez, Mon. Not. Roy. Astron. Soc. 404, 885 (2010). 



15 



